Evidence for the E 2u model of the superconducting order parameter 
in UPt 3 from Josephson interferometry 
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We present data on the modulation of the critical current with applied magnetic field in UPt3-Cu- 
Pb Josephson junctions and SQUIDs. The junctions were fabricated on polished surfaces of UPt3 
single crystals. The shape of the resulting diffraction patterns provides phase sensitive information 
on the superconducting order parameter. Our corner junction data show asymmetric patterns with 
respect to magnetic field, indicating a complex order parameter, and both our junction and SQUID 
measurements point to a phase shift of 7r, supporting the E^u representation of the order parameter. 



INTRODUCTION 

More than two decades after its discovery PQ, the mech- 
anism of superconductivity in UPt3 is still unknown. Al- 
though UPt3 has many unusual properties which suggest 
that the order parameter is unconventional, the pair- 
ing symmetry has not been unambiguously determined. 
Perhaps most unusual is that it exhibits a characteris- 
tic double peak in specific heat [2 , indicating two dis- 
tinct superconducting phases, with an initial transition 
at T c+ w 5b0mK and a second transition T c _ w 500mK. 
Subsequent measurements revealed a third phase at high 
magnetic field [3 j. Transport measurements show power 
law dependencies at low temperatures, revealing the pres- 
ence of nodes in the superconducting gap[H [5]. Muon 
spin resonance showed signs of spontaneous magnetiza- 
tion, and thus time-reversal symmetry breaking (TRSB) 
in the low temperature phase [6], but this result has not 
been reproduced in later measurements [7 . NMR stud- 
ies of the Knight shift support a triplet pairing mech- 
anism or possibly a singlet state with strong spin-orbit 
scattering [8]. 

Numerous models for the pairing symmetry have been 
put forward to explain this complicated behavior, but the 
two best candidates are the singlet-state E\ Q and triplet- 
state E2 U representations of the order parameter [9j HO] . 
Both models feature a real order parameter in the high 
temperature A-phase and a complex order parameter be- 
low the second transition in the B-phase. Various efforts 
have been made to distinguish between these two theories 
on the basis of experiment [TTJ [12j [T3l [14] , in particular re- 
lying on details of transport properties, but the relatively 
subtle differences in node structure and gap magnitude 
have proven difficult to resolve. Perhaps the clearest dif- 
ference between these models is the periodicity of phase 
winding in the order parameter, as seen in Table 1. A 
rotation of 90° about the c-axis causes a phase shift of 
7r in the E\ g model, but a phase shift of 2tt in the E2 U 
model. In this paper, we propose to distinguish between 
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TABLE I: (color online) Graphical depictions of the E\ g and 
E2u models of the order parameter for UPt3. Columns de- 
note the high and low temperature (A and B, respectively) 
superconducting phases, and rows denote the two theoretical 
models. 



these two models by detecting this difference in phase 
with Josephson interferometry. 

Josephson interferometry, used successfully to charac- 
terize the cuprates as d-wave and Sr2Ru04 as complex 
p-wavep~5j [T6j [17], remains the most definitive phase- 
sensitive test of the order parameter of unconventional 
superconductors. In this technique, a superconducting 
weak link is created between two superconductors - in 
our case, a single crystal of UPt3 and a film of the con- 
ventional superconductor Pb. Applying a magnetic field 
to this junction perpendicular to the current flow creates 
a phase gradient along the junction that alters the lo- 
cal supercurrent density. Any intrinsic phase differences 
arising from the order parameter symmetry will also af- 
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feet the current density. The critical current (I c ) can be 
given as a function of external flux (& ex t) and intrinsic 
phase difference (5) by the following: 



sin (jr^ext/^Q + S/2) 



(1) 



In the case of uniform s-wave superconductors, this re- 
sults in the conventional Fraunhofer diffraction pattern 
shape for plots of critical current vs. applied magnetic 
field. Even in a superconductor with an anisotropic order 
parameter, as long as the junction is on a single flat crys- 
tal face, the pattern will look Fraunhofer. This is because 
tunneling probability falls off exponentially with barrier 
thickness and so the Josephson current effectively probes 
a single k-space direction. If, however, the junction wraps 
around the corner of a superconductor with a sign change 
in the order parameter between the two tunneling di- 
rections, part of the junction will probe each direction, 
and the phase shift will cause a distinctive change in the 
diffraction pattern. The predicted patterns for each of 
the order parameter models are given in Figure [I] 

EXPERIMENT 

The UPt3 crystals were grown in an electron-beam 
floating zone furnace, and have a measured residual re- 
sistivity ratio (RRR) of greater than 900, in some cases 
as high as 1100, indicating their exceptional purity. We 
polished the surfaces with 0.3/im diamond lapping films, 
and then glued them to a glass substrate with Pyralin® 
polyimide coating. After masking with a dry photoresist, 
the were ion milled and 150nm of Cu was evaporated as 
a normal metal barrier, followed by 800nm of Pb as the 
superconducting counter-electrode. Junction dimensions 
were typically 50xl00/im. Junctions on UPt3 are known 
to be susceptible to magnetic flux trapping [18 so the 
samples were cooled in a Kelvinox® dilution refrigerator 
with Cryoperm® and lead cans to provide the necessary 
magnetic shielding. The junction voltages were excep- 
tionally small (« nV), and so were measured with a 
superconducting quantum interference device (SQUID) 
potentiometer circuit, in which an inductively-coupled 
SQUID detected the current flowing through a known 
resistor in parallel with the junction. 

All the measurements reported here were carried out 
in the low temperature B-phase of UPt3. The junctions 
exhibited nearly ideal resistively-shunted junction (RSJ) 
behavior, as well as showing Shapiro steps when an ac 
modulation was applied. Junctions fabricated on a sin- 
gle crystal face displayed diffraction patterns that were 
nearly Fraunhofer and symmetric around zero field, indi- 
cating uniform phase and no trapped vortices. Examples 
of these measurements, as well as a sample photo, can 
be seen in Figure [2] It is worth mentioning that even 
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FIG. 1: (color online) Planar representations of the order 
parameter laid on top of a schematic of a corner junction, 
with the corresponding diffraction pattern placed alongside, 
(a) An s-wave order parameter produces the classic Fraunhofer 
pattern. (b)The E\ g B-phase produces an asymmetric double 
peak. (c)The E^ u B-phase produces a symmetric double peak. 



though the B-phase of UPt3 is expected to be chiral and 
exhibit TRSB, similar to Sr 2 Ru04, we saw none of the 
evidence for chiral domains in UPt3 that were seen in 
Sr2Ru04[17 , such as hysteresis or switching noise. 

We also measured junctions fabricated so that they 
straddled the corner between the a-and b-axes. These 
corner junctions behaved quite differently than the edge 
junctions. In all cases the diffraction patterns they pro- 
duced exhibited features that were asymmetric with re- 
spect to field polarity. The asymmetry was not caused 
by the self- field effect of the current through the junc- 
tion, because the pattern was symmetric with respect to 
the direction of current flow. Changing temperatures af- 
fected the magnitude of the critical current, as expected, 
but left the shape of the patterns unchanged. This asym- 
metry with respect to field polarity is a characteristic sign 
of a complex order parameter and TRSB, which supports 
both models' predictions of a chiral state in the B-phase. 

Aside from the ever-present asymmetry, the corner 
junction patterns varied between different junctions or 
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FIG. 2: (color online) (a) A photograph of one of our samples 
- the crystal is the large block in the upper middle with four 
edge junctions evaporated on its surface. The lines of Pb 
film and indium pads are also visible. (b)A characteristic 
IV plot for a junction, showing classic RSJ behavior. (c)A 
diffraction pattern for one of our edge junctions, exhibiting a 
nearly Fraunhofer shape, (d) Shapiro steps from a junction, 
confirming Josephson behavior. 



thermal cycles of the same junction, with three or four 
qualitatively similar patterns recurring. The variations 
in the patterns between different junctions could be due 
to surface roughness at the corner. It is unlikely that the 
two polished surfaces make a perfect corner, and faceting 
at the region where they meet is probable. These facets 
can cause distortions in diffraction patterns, particularly 
for non-s-wave superconductors, but even in this case an 
asymmetric variation is a signature of a complex compo- 
nent of the order parameter [T9 l [20] . 

The changes in the diffraction patterns between ther- 
mal cycles require a dynamic mechanism to explain. The 
most obvious candidate is vortex trapping near the junc- 
tions. Though our magnetic shielding (H res « 10 _4 G) 
combined with very slow cooling cycles seems to be suf- 
ficient to prevent vortex entry into the edge junctions, it 
is possible that the corners of our samples are damaged 
enough to provide more pinning sites. With this in mind, 
we have tried modeling corner junctions combining an in- 
trinsic phase shift with trapped vortices. We tested phase 
shifts corresponding to the three candidate symmetries: 
(s-wave), tt/2 (E\ g ), and tt (E2 U )- We modeled a vor- 
tex as a Gaussian contribution to the flux through the 
junction with integrated flux = <£ /2 and width equal to 
3-5% of the junction width. We then tried to reproduce 
the types of patterns we measured most often. We do not 
claim to have modeled the junctions exactly, but the best 
qualitative fit to the observed patterns was given by an 
intrinsic phase shift of tt with vortices close to the cen- 



ter of the junction. Comparisons of a few representative 
diffraction patterns to simulations are given in Figure [3] 
In an effort to avoid the complications caused by the 
material properties of the corners, we fabricated two 
junctions, one on either side of the corner, incorporat- 
ing the crystal into a SQUID geometry. In this case, an 
intrinsic phase difference in the crystal will show up as 
a shift in the peak of the critical current modulation, as 
given by: 



I c {$ext) = 2J 



(2) 



Unlike a single junction, the periodic modulation of 
a SQUID does not provide a central peak to reveal the 
zero point of magnetic flux, so extra care is required to 
rule out extrinsic shifts in the pattern. If the two junc- 
tions are not identical, unequal current flow between the 
branches of the loop will couple field into the SQUID, 
shifting the pattern. To account for this, we bias the 



Data 



Simulation 



(a) 


t 


0.8 
0.6 








5 °- 4 

0.2 
0.0 


aA/UAa 



-250 250 

Field (mG) 




-500 -250 250 500 

Field (mG) 



FIG. 3: (color online) Comparisons of corner junction diffrac- 
tion patterns with simulations. All simulations shown as- 
sume the Em representation, (a) Corner junction with no 
vortex trapping, (b) Junction with a vortex trapped within 
5% of junction width of the corner, (c) Junction with a vortex 
trapped within 15% of junction width of the corner. 
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FIG. 4: (color online) (a)Two SQUID (one junction on either 
side of a corner) modulation curves taken at different bias 
currents - arrows denote the peaks and highlight the shift 
in position caused by asymmetric current flow. The arrow 
locations for these and other curves correspond to data points 
in plot (b). (b) Extrapolations to zero bias current for seven 
thermal cycles of the corner SQUID. The lines cluster around 
a phase shift of 0.5<3>o. 



SQUID at various current levels as shown in Figure [4^i, 
noting peak location at each current, and then extrapo- 
late the peak location to zero bias current [21] . In order to 
account for residual external field or trapped flux, we per- 
formed several thermal cycles, with the results plotted in 
Figure ^p. The results cluster around 3>o/2, correspond- 
ing to a phase shift of n, which agrees with our corner 
junction results and also supports the E2 U model. 



CONCLUSIONS 

In summary, we have fabricated Josephson junctions 
on high quality single crystals of UPt3 and used them to 
perform phase sensitive measurements of the supercon- 
ducting order parameter. There is strong evidence for 
a complex component of the order parameter from the 
asymmetry with respect to field polarity in the diffraction 
patterns, but we found no sign of chiral domains in our 
edge junctions. By comparing our corner junction results 
with simulations involving vortex trapping at the corner 
and from shifts in SQUID modulation curves, we find ev- 
idence for an intrinsic phase shift of tt for a 90° rotation, 
in agreement with the Em representation of the order pa- 
rameter. We are continuing further measurements with 
different surface treatments, including as-grown crystal 
faces, in order to reduce the effect of vortices, as well as 
studying the crossover between the A and B phases. 
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